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In this paper, we compare four different methods of fault simulation 
in terms of their handling of arbitrary numbers of logic values, 
modeling levels, and detailed timing. The methods considered are 
parallel, deductive, multilist, and concurrent simulation methods. 
Since some of the methods, in their current forms, are unable to 
handle all the problems under consideration, we have proposed 
extensions to the methods wherever necessary before making the 
comparisons. While all the methods considered are capable of solving 
the problems with the same degree of accuracy, the concurrent sim- 
ulation method appears to be the simplest and most flexible. 

I. INTRODUCTION 

Different techniques for the efficient simulation of faults in digital 
circuits have been published. Among these, the best known are parallel 
simulation, 1 " 3 deductive simulation, 4 and concurrent simulation. 5 A few 
papers analyzing some aspects of these methods have also been pub- 
lished. 6 " 9 

This paper and two others 1011 comprise a series attempting a com- 
prehensive analysis of fault simulation methods. It is hoped that they 
will provide a basis for the selection of fault-simulation methods to 
satisfy specific requirements. 

In this paper, we consider three aspects of circuit modeling and their 
effects on the fault-simulation method used. First, we consider the 
number of logic values needed to accurately model logic devices and 
its impact on the simulation method. Next, the effectiveness of the 
different methods for simulating at different levels (e.g., gate level, 
functional level, subsystem level, etc.) is considered. Finally, we discuss 
the modeling of timing effects, such as rise and fall times and high- 
frequency rejection. 

Our study covers four methods of fault simulation: parallel, deduc- 
tive, multilist, and concurrent. In their current forms, some of the 
methods are not capable of handling all the problems we consider. 
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Therefore, we have attempted to extend the existing methods, wher- 
ever necessary, before making the comparisons between methods. 
Before proceeding to the analysis of the methods, we present a brief 
description of each method. 

Historically, parallel simulation was the first method that simulated 
a number of faults simultaneously. 1 This method, which is perhaps the 
most widely used, takes advantage of word-oriented operations in the 
host computer and packs together several faulty circuit values into one 
or more computer words. Although this method is quite efficient, 
multiple passes are required for simulating large numbers of faults. 

Deductive simulation attempts to eliminate the need for multiple 
passes by computing normal signal values in the circuit and deducing 
the faulty values by manipulating lists of faults. 4 Associated with each 
signal is a fault list, which is a set of faults, any one of which will cause 
the signal value to be different from the normal value. The effects of 
faults are propagated through the circuit by an algebra of sets. 

The multilist method associates two or more lists of faults with each 
signal. 1012 Conceptually, the number of lists associated with a signal is 
equal to the number of logic values simulated. Thus, for two logic 
values, there will be a O-list and a 1-list associated with each signal, 
the former being the set of faults in whose presence (individually) the 
signal will have the value 0, and the latter those that result in a value 
of 1. Set algebra is necessary for manipulating these lists also. However, 
unlike the deductive method, the equations for computing the output 
lists of a device from its input lists are dependent only on the function 
performed by the device and not on the signal values. 

In concurrent simulation, any fault that causes the inputs, outputs, 
or internal state of a device to be different from their normal values is 
represented conceptually by a copy of the device. During simulation, 
if the inputs, outputs, and state of a faulty copy become identical to 
those of the fault-free copy, the faulty copy is deleted. Thus, faulty 
copies are created and deleted during simulation. The evaluation of 
faulty copies is essentially the same as fault-free copies, and no set 
algebra is involved. Concurrent simulation can also handle a large 
number of faults simultaneously. 

It is interesting to note that all the above methods, except parallel 
simulation, use some form of data compression for storing faulty signal 
values. On the other hand, parallel simulation attempts to compute 
simultaneously the fault-free signal value and a number of faulty signal 
values associated with each lead in the circuit. 

II. NUMBER OF LOGIC VALUES 

Three-valued logic systems have been widely used for analyzing 
essentially binary systems. 1314 Three logic values are also used in logic 
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simulation, where and 1 represent the two discrete values being 
modeled and a third value, u, denotes that a particular value is 

unknown. 

Recently, tri-state busing has become a widespread technique used 
in many lsi designs. Difficulties in modeling effects associated with 
CMOS technology have been reported. 15 One effect is the memory 
associated with a disabled bus. That is, the disabled bus remembers 
the previous logic value on the bus. A solution consists of adding 
special circuitry to regular gates, making possible the use of a simulator 
with only three logic values. 15 An alternate solution is the addition of 
three more logic values, namely z , z u and z, for representing the states 
of disabled buses, with previous value equal to 0, 1, and unknown, 
respectively. 16 Transistor-transistor logic (ttl) tri-state technology 
requires the addition of only one logic value, 2. 16 

Bus contention, another typical, potentially destructive tri-state 
effect, cannot be modeled by added circuitry. A solution consists of 
adding one more logic value representing a conflict state, a, as shown 
in the following example. 

Consider a driver inverter and a bus configuration in ttl tri-state 
technology (Fig. 1). When line e is enabled, the gate operates as an 
inverter, when e is disabled the output of the gate is in a high- 
impedance state. When used in a bus configuration, two enabled 
inverters create a conflict (bus contention), if they are in opposite 
states. The set of logic values {0, 1, u, a, z) is sufficient to model these 
effects, since tri-state devices in ttl technology do not have the 
memory property mentioned above. 

Table I shows how the bus configuration of Fig. 1 can be simulated 
using the above set of five logic values. Since the bus will be connected 
to the output of drivers, which can produce four out of the five logic 
values, only four logic values are used for modeling the bus. 




(b) 
Fig. 1— (a) TTL Driver-inverter, (b) Bus configuration. 
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Table I — (a) Tristate inverter output 
(b) State of tristate bus 
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If an ordinary gate could be connected directly to a bus, the model 
should allow five logic values for the gate inputs, but requires only 
three logic values for its output. Table II shows the behavior of such 
an AND gate with inputs x and y, and output t. 

The use of larger sets of logic values, though necessary to correctly 
model modern technology, has a serious impact on the method of 
simulation used. The following sections deal with this problem. 

2.1 Parallel simulation 

When using a switching algebra (i.e., two logic values) parallel 
simulation can be implemented by associating one computer word 
with each line in the circuit. One bit of this word represents the signal 
value on a line in the fault-free circuit and the remaining bits represent 
values on the same line in the presence of different single faults. 

Table II — AND gate with five input 
logic values 
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Table III — Coding for three logic 
values 
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When a three-valued system is used, each of the circuits simulated 
in parallel must be coded using two binary digits. A commonly used 
method consists of associating two words with each line a, namely the 
0-word, a , and the 1-word, a 1 . 17 The coding used is shown in Table III, 
where the subscript i refers to the ith bit of each word. 

Examples of its use are shown in Fig. 2. Here, and elsewhere in this 
paper, lower-case roman letters are used to denote leads and Greek 
letters represent words. For the gates of Fig. 2, we have 

y° = a + P° 

6° = a°.^° 
5 1 = a 1 + P 1 

where • and + represent the bitwise AND and OR operations on 
complete words. 

This method can be extended for any number of logic values. For 
instance, consider the AND gate of Fig. 3, using the set of logic values 
{0, 1, u, a, z). The binary coding scheme requires three computer 
words for each line, and three codes (out of eight) are not used. For 



■<:■ 



AND 






(3° 



< t.0 




OR 



Fig. 2 — Use of coding to represent signal values on gates. 
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Fig. 3 — AND gate representation for five logic values. 

any choice of code, it is possible to calculate the gate output from 
switching expressions of the following form: 

Y° = f(a°, a\ a 2 , 0>, fi\ ?) 

y 2 = h(a ,a\a 2 ,(?,fl\fl 2 ). 

The original set of logic values and operations do not constitute a 
Boolean algebra. The coding scheme establishes a mapping of non- 
Boolean functions into switching operations that can be applied on full 
computer words, thus, allowing parallel processing. 

By using a coding of n — 1 variables to represent n logic values, it is 
possible to obtain simpler equations for computing the outputs of 
gates. For example, consider the gate of Fig. 3 and a coding using four 
words a , a 1 , <x a , and of to represent five logic values. The code is such 
that a J i = 1, j = 0, 1, a, or z, iff a, = j. All the variables will be zero, if 
and only if a, = u. With this coding, the following equations are 
obtained for the gate of Fig. 3: 

y° = 

This type of coding can be used for any number of logic values. 

2.2 Multilist simulation 

It has been shown that for a three-valued logic system, three lists 
X°, X 1 , and X u can be associated with each line *. 1012 Each list, X 1 , 
(i = 0, 1, u) represents the faults which cause line x to have the value 
i. For each line x, all the lists X 1 are disjoint and any list is the 
complement of the union of the other two (i.e., the union of the three 
lists is the set of all faults being simulated). 

For the gates of Fig. 2, we have 
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C° = A UB Q 

D l =A 1 L)B 1 
D° = A°n B° 

X 1 = A 



c = c° u c 1 



D" = DP U D l 



yO _ Al X — X U X , 

where ~, U, and n, are set complement, union, and intersection, 
respectively. 

When five logic values are used, we need five lists; for instance, A , 
A\ A", A", and A z are associated with line a. 

For the AND gate of Fig. 3, we have 

C° - A U B° 
C 1 =A l D B l 
C a ={ } 
C*={ } 



C U =(C°UC'UC'U C a ) =i UB°UA 1 flB 1 . 
For the inverter of Fig. 1, we have 
B° = E 1 n A 1 
B 1 = E i D A 
B* = E° 



B u = E" U E a U E z U (E l D (A U A 1 )) 
B"={ ), 
and for the bus configuration of Fig. 1 

A° 3 = (B? n B° 2 ) u (S? n Bl) u (S^ n B\) 
Al = (B\ n B2 1 ) u (fi} n 51) u {B\ n Bf) 

AS = Bt U B^ 
Ai = Bt n BI 

A? = (B° 2 n Bl) u (Bl n B?). 

This method can be generalized to any gate type and any number of 
logic values as follows: Let us assume that we wish to simulate a 
function f(x h x 2 , • • • , x n ), where each input and the output may 
assume any one of k values, denoted by 1, 2, • • • , k, and that the 
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function is defined by a table which specifies the values of f for all 
combinations of values of x»-. 

(i) We associate a variable x\ with each variable Xi, such that 
x\ '■= 1 if and only if x, = j, 1 < y < &. Similarly, we associate A variables 
Z 7 ' with /. 

(ii) For each i, 1 < i < &, we obtain an expression 

where P, are products of literals xf, representing all combinations of 
values for which f = i. For example, if the table has an entry 

Xi = 1, & — 0, X3 — z, f = 1, 

the expression for f l will contain the term 

X1X2X3. 

(Hi) Replace all lower-case letters in the equation for f by the 
corresponding upper-case letters, representing lists, and retain the 
superscripts and subscripts. Replace products by intersection and sums 
by union. 

2.3 Deductive simulation 

Deductive simulation is well defined for two logic values, and is also 
applicable to three logic values with some loss of information. 4 Specif- 
ically, if the signal value in the normal circuit is known, (i.e., or 1), 
but the value in the presence of a fault a is unknown (denoted by u), 
the fault a is included in the fault list as a star fault; 418 that is, it is 
unknown whether the particular signal value in the presence of the 
fault a will be different from the fault-free value. It was shown in Refs. 
10 and 12 that there are cases where the circuit value in the normal 
circuit may be unknown, but the value in the presence of a fault may 
be known. Since the deductive method cannot represent this case, the 
results obtained may be less accurate than with other methods. 10,12 

A modification of the deductive method that leads to accurate three- 
valued fault simulation was presented in Ref. 10. It uses the coding of 
Table III for representing each signal value by a pair of binary 
variables. A pair of equations can then be derived, as in Section 2.1, 
for computing the coded outputs for each gate type. These equations 
can be viewed as defining a transformation of the original circuit with 
three signal values into two circuits that will have only binary signals. 
These two circuits can be simulated using the two-valued deductive 
method. The fault-free and faulty signal values on any lead in the 
original circuit can be determined from the signal values and fault lists 
associated with the corresponding pair of leads in the transformed 
circuits. 
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Fig. 4 — Tristate bus driver. 

The same approach can be used for performing deductive fault 
simulation with any number of logic values. If k logic values are to be 
simulated, flog2£l binary variables will be used to represent them, 
where \x\ denotes the smallest integer greater than or equal to x. The 
equations for the coded outputs of different gate types can be derived 
from their truth tables, and used in deductive simulation. 

As as example, consider the bus driver of Fig. 4 to be simulated with 
four logic values, namely, 0, 1, z (high impedance) and u (unknown). 
The behavior of the device is specified in Table IV. 

Using the coding of Table V, we shall represent the signals a, e, and 
b of the bus driver by ao and ai, e and e h and b and 6i. The output 
equations 6 and b\ can be derived from Tables IV and V. 

b = e -ei + Oo-eo-ei 

bi = eo-ei + Oi-Co'Ci 

For any combination of input values and fault lists, the output values 
and fault lists can be computed as in Ref. 19. 

Denoting the fault list associated with each variable by the corre- 
sponding upper-case letter, let the input values and fault lists for the 
circuit of Fig. 4 be as follows: 



ao = 0; 


Ao = {1, 3} 


ai-1; 


A, = {3} 


e = 0; 


E = {2, 4} 


e, = 1; 


Bt = {4, 5} 
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Table V — Coding for tristate devices 
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Let us assume that all the faults being considered are external to the 
device, and we wish to propagate the effects of the faults through the 
device. The input conditions are: a = 1, e = 1. Since the fault 1 is 
contained only in the fault list Ao, it will cause ao to become 1, and 
therefore a to become z. On the other hand, fault 3 is contained both 
in Ao and A u and will cause both ao and a i to be inverted; the value of 
a in the presence of fault 3 will be 0. Similarly, fault 2 will result in 
e = z, fault 4 in e = 0, and fault 5 in e = u. 

For the above set of values, the output values and fault lists can be 
computed using the equations for 6o and 6i and the method presented 
in Ref. 19 as follows: 

6 = 

&! = 1 

Bo = (Ao n W n ~Ei) U (E n Ex) = {1, 3, 4} 



Bx = (Ai U Eo U Ei) H (Eo D Ei) = {2, 3, 5} . 

Denoting the value of b in the presence of fault a by b(a), we can 
obtain the following faulty values from the values b and bi and fault 
lists Bo and B\. 

b(l) = z; 6(2) = u; b(3) = 0; 6(4) = z; 6(5) = u. 

These can be verified by computing the output for each faulty com- 
bination of inputs using Table IV. 

The modified deductive method discussed above does not lose any 
information about the normal and faulty circuits and is as accurate as 
any of the other methods. It requires only |"log2/il lists compared to the 
n lists needed for multilist simulation. However, fault list computations 
depend on signal values and may be more complex than in the multilist 
method. 

2.4 Concurrent simulation 

There is no limitation on the number of logic values in this simula- 
tion method since faulty and fault-free circuits are treated independ- 
ently. As long as the primitive elements of the circuit are well defined, 
the evaluation of faulty circuits presents no difficulty. 
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2.5 Summary of results 

The results of Section II are represented in Table VI. 

Deductive simulation with three logic values (indicated by * in 
Table VI) requires the introduction of the concept of star faults. 
Deductive fault simulation for more than three logic values (indicated 
by f in Table VI) could be defined by using a transformed circuit as 
proposed in Section 2.3. However, the complexity of such a procedure 
does not seem to justify its use. 

From the point of view of simulating more than three logic values, 
concurrent simulation represents the simplest, most flexible simulation 
method. 

III. MODELING LEVELS 

Three levels of modeling and their effects on the simulation method 
used will be considered: gates, higher-level primitives, and user-defined 
functions. 

3. 1 Gate-level simulation 

All the fault-simulation methods presented here were initially de- 
veloped for simulating circuits modeled at the gate level. Therefore, 
none of the methods presents any problem, provided only two (or 
three) logic values are to be simulated. The differences due to the 
number of logic values needed have already been discussed in Section 
II, and the effects of detailed timing analysis are discussed in Section 
IV. 

3.2 High-level primitives 

It is often convenient to model devices such as flip-flops, multiplex- 
ors, counters, and shift registers as high-level primitives rather than as 
interconnections of gates. For purposes of simulation, such devices 
may be described by tables, Boolean equations, or algorithms. The 



Table VI — Summary of Results: logic values 
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simulation with three logic values. 

fault simulation for more than three logic values. 
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type of representation that is most convenient to use will usually 
depend on the simulation method. 

3.2. 1 Parallel simulation 

Several solutions are possible. The input values for individual faults 
can be determined from the input word(s), and the outputs of the high- 
level primitive can be evaluated for each case. The output values must 
then be packed so that the parallel simulation method may be used 
elsewhere. While this approach may be satisfactory for predominantly 
gate-level circuits which also contain a few high-level primitives, the 
overhead associated with converting to single-fault simulation and 
back to parallel simulation may not be acceptable. 

When only two logic values are involved, the primitive can be 
represented by Boolean (switching) expressions. The operators in 
these expressions can be treated exactly like gates in parallel simulation 
of gate-level circuits. When more than two logic values are simulated, 
the description of the primitive may be in the form of tables. Using a 
coding of the type described in Section II, switching algebraic expres- 
sions for the coded output words (i.e., the 0-word, 1-word, etc.) can be 
obtained in terms of the coded words associated with the inputs and 
state variables. These equations can then be used to compute the 
coded output words. 

3.2.2 Multilist simulation 

The function realized by a high-level primitive can be represented 
by tables. From these tables, equations for the output lists in terms of 
lists associated with inputs and state variables can be obtained as 
discussed in Section 2.2 and used for simulation. 

3.2.3 Deductive simulation 

As in the case of parallel simulation, one approach is to simulate 
each high-level primitive for one fault at a time and use the results to 
construct output fault lists for use outside the primitive. Alternatively, 
outputs and the next state values of state variables may be represented 
by Boolean equations, which are used for fault-list computations in 
the same manner as gates. When more than two logic values are to be 
simulated, a binary coding can be used as discussed in Section 2.3, and 
equations for each coded bit can be used for fault-list computations. 

Another possibility is to use tables that specify the fault-list com- 
putations for every combination of input values. 4 These tables can be 
constructed from the tables specifying the primitive, as shown by the 
example of Fig. 5. 

The behavior of the function is represented in Table VII, where p h 
qi are the initial states of the flip-flops and p 2 , 92 are the next states. 
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Fig. 5— NAND sr latch. 
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-Behavior of SR latch 
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The table does not include 00 as an initial state because it cannot be 
produced directly. 

The fault-list propagation is summarized in Table VIII and does not 
include local faults. Pi, Qi, R, S, P>, and Q 2 are the fault lists associated 
withpi, <7i, r, s, p 2 , and q 2 . The star faults in the table are to be added 
to both the fault lists, P 2 and Q 2 . We shall demonstrate the procedure 
used for deriving Table VIII, by showing how line 1 of the table was 
obtained. 

Consider line 1 of Table VII. To get a change inp 2 , we need a change 
into lines 3, 4, 7, or 11 and for a change in q 2 , we need a change into 
lines 2, 6, 8, or 10, which produces: 

P 2 = P i QiRS U PiQiRS U P^RS U PiQiRS U P,Qi(4>) 

= RS\J PiR 
Q 2 = P&RS U P^RS U PiQiRS U PiQ^S U PiQi(<S>) 

= RS u Q,S, 

where juxtaposition represents intersection. Since a change to line 12 
is needed to cause p 2 and q 2 to become unknown, the star faults are 
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Table VIII— Fault-list equations for SR latch 
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given by 



PiQjlS 



We have used PiQi, which corresponds to the initial state 00 in the 
faulty circuit, as a don't-care state (O) to simplify the expressions. 

3.2.4 Concurrent simulation 

In concurrent simulation, the same method is used to evaluate fault- 
free and faulty circuit signal values. Therefore, no transformations of 
representation are necessary, and any representation that leads to 
efficient simulation may be chosen. 

3.3 User-defined functions 

Our discussion of Section 3.2 also applies to user-defined functions. 
The main difference is that the tables or equations used for represent- 
ing the functions must be generated from descriptions in a high-level 
language such as the function definition language in LAMP. 20 

A typical construct in such a language is the cause-effect statement. 
Such statements can be nested to many levels. The techniques dis- 
cussed in Section 3.2 can be used for simulating user-defined functions 
using the parallel, multilist, or deductive method by first replacing 
cause-effect statements by equivalent equations. For example, the 
statement 

if x then z = a else z = b 



can be replaced by 



z = a-x + b'X + a-b. 



The redundant term a-b has been introduced to produce the correct 
result z = 1 for the case a=6 = landx = u. 21 Otherwise, the 
pessimistic result z = u will be produced for this case. 
Concurrent simulation does not require the transformation of cause- 
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effect statements into equations. For each fault and each combination 
of inputs and state, only those computations enabled by the conditions 
need be performed. The operations in a function definition need not 
be restricted to logical operations. Therefore, it is not necessary to 
generate Boolean equations corresponding to arithmetic operations, as 
would be necessary in the other methods considered. Thus, it appears 
that the concurrent method would allow simulation of functions de- 
fined at a higher level than is possible with the other methods. 

3.4 Summary of results 

The results of this section are summarized in Table IX. The concur- 
rent method is clearly superior in its ability to simulate different levels 
of models. 

IV. TIMING 

In this section, we study different effects related to timing, and their 
impact on the four simulation methods under consideration. We shall 
consider the effects of different rise and fall times associated with 
signal changes, suppression of short pulses to model inertial delays, 
and the simulation of faults which affect the magnitude of delays 
associated with devices. We shall restrict our discussion to logic 
simulation with two and three logic values. 

4. 1 Rise and fall times 

The delays associated with to 1 and 1 to transitions of a signal, 
called here the rise and fall times, are not necessarily equal. 22,23 All the 
methods of fault simulation under discussion simulate a number of 
signals simultaneously, some of which may be rising and some falling. 
To simulate this effect accurately, a mechanism is necessary for 
allowing rising and falling signals to change at different times. 

Let fobea time before any change occurs on the line under consid- 
eration. Due to differences in the rise and fall times, signal changes 
may occur on the line at times t\ and t 2 , where t < t\ < t 2 . Thus, at 
time t 2 , all signal changes associated with the particular event would 
have occurred. The effect of different rise and fall times can be 

Table IX — Summary of results: modeling levels 





Parallel 


Multilist 


Deductive 


Concurrent 


Gate level 


1 


1 


1 


1 


Higher level 
primitives 


2 


2 


2 


1 


User defined 
functions 


2 


2 


2 


1 



Note: 1 = No transformations required. 

2 = Transformation into equations required. 
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simulated accurately by computing the values of the signals at time t\ 
based on the values at t and t 2 , namely, the initial and final values for 
the particular set of transitions. 

In the following sections, we shall consider four simulation methods 
and examine the results produced by their different models at three 
points in time, namely, before ti(to), between ti and t 2 {t\), and after 
fe(fe). 

4.1.1 Parallel simulation 

Let £o, £1, and £> be the words associated with a line x at times to, t u 
and t 2 , in two-valued parallel simulation. In three-valued simulation, 
two words denoted by superscripts and 1 will be associated with the 
line for each of the above times, and the coding of Table III will be 
used. 
Case 1: Rise time < fall time. We have 

£i = & + £> 

for two-valued simulation, and 

e - a + g 

for three-valued simulation where + and • represent bitwise OR and 

AND performed on full words. 

Case 2: Fall time < rise time. We have 

for two-valued simulation, and 

/rO _ t° j. t° 
«1 - SO + S2 

for three-valued simulation. 

The preceding formulas can be verified by checking all nine possible 
transitions between the set {0, 1, u) and itself. 

4.1.2 Three-list method 

Let Xo, X[, and X 2 be the i-lists at times to, fe, and t 2 , defined above, 
for i = 0, 1, u. Using the same arguments as in Section 4.1.1, we obtain 
the lists for time t\ as given below. 
Case 1: Rise time < fall time. We have 

X\ = Xo u x 2 
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Case 2: Fall time < rise time. We have 

X\ = .Xo n x.2 
x\-x%uxl. 



In both cases, X\ = (X? U X\). 

4.1.3 Deductive simulation 

Deductive simulation with different rise and fall times has been 
discussed by Kjelkerud and Thessen. 24 Here we present an alternate 
method. 

Let the times to, t\, and tt be as defined earlier and let Xi and Xi 
represent the signal values and fault lists at those times, i = 0, 1, 2. If 
the rise time is less than the fall time, all to 1 transitions will occur 
at t\. Therefore, we have 

X\ = Xq + Xo . 

Similarly, if fall time < rise time, 1 to transitions will occur at t h and 
Xi = 1 if and only if it remains at 1 throughout the transitions. 
Therefore, for this case 

X\ = ATo • J^2 • 

The fault lists X\ at time t\ for different signal changes in the fault-free 
circuit and different relative values of rise and fall times can be 
determined from these equations. They are summarized in Table X. 

4.1.4 Concurrent simulation 

In fact, concurrent simulation is a trivial case, because fault-free and 
and faulty circuits are simulated independently. Rising and falling 
edges will still occur in distinct event waves, but the treatment of these 
events is individual. 

4.2 High-frequency rejection 

High-frequency rejection consists of eliminating short pulses for 
modeling the effect of inertial delays. We consider events occurring at 

Table X — Fault-list equations for handling different rise and fall times 



Fault-Free Circuit 


Rise Time < Fall Time 


Fall Time < Rise Time 


Rising edge 

Xo — 0; X2 = 1 


x, = Xo n x- 2 


x,=0_ 

x,=x nx 2 


Falling edge 
Xo = 1; X2 = 


•r, = 1 _ 
A i — Ao ii A 2 


X,_z=0 

x, = Xo n x 2 


Constant one 
*o ■ 1; X2 m 1 


x, = l 

x, = Xo n x 2 


Xi«l 

Xi = Xo U X 2 


Constant zero 
Xo = 0; Xi ■ 


xi -0 

X, = Xo u x 2 


x, = 

A i = An 11 A2 
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times to, t h and t 2 , where t < fa < t 2 and the logic values at these times. 
If h - fa is less than the magnitude of the inertial delay, then the 
change at fa must be rejected to suppress short pulses and the value 
between U and h, x u will be replaced by a corrected logic value, Xm. If 
xo, Xi, and x 2 are the computed signal values at to, fa, and h respectively, 
and if h - fa is less than the inertial delay, the corrected signal value 
at time fa is given by: 

x-i n = X0X1 + x\X 2 + XoX 2 . 

The method for performing high-frequency rejection can be derived 

from this equation. 

In case there are more than two events within the range of the 
inertial delay, the treatment elaborated above must be repeated for 
each pair of events within that range. For example, consider three 
events occurring at times fa, fa, and fa. The following triples will be 
considered: (to, fa, t 2 ), (t , t u fe), (*i, t 2 , t 3 ), which represent three pairs 
of events. 

4.2. 1 Parallel simulation 

For two-valued parallel simulation, the word & has to be replaced 
by 

£i« = loli + && + &6 

and for three-valued parallel simulation, £? and |! are replaced by 

& = loii + m + 10I2 1 

The coding denned in Table III was used to obtain the above equations 
for three-valued parallel simulation. 

4.2.2 Three-list methods 

Using the equation for x in given above, we obtain the following fault- 
list equations: 

x° ln - (X° n x?) u (X? n x° 2 ) u (X° n x° 2 ) 
x\ n = (XI n xl) u (X! n x\) u {Xl n x\) 



xt n = (X? u XI). 

4.2.3 Deductive simulation 

The deductive fault list X ln can be obtained from the equation for 
the new signal value Xm, the values of Xo, Xi, and x 2 , and the associated 
fault lists. The fault list X ln can be computed in the same manner as 
fault propagation through functional blocks. 19 In fact, high-frequency 
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rejection may be thought of as being performed by a filter whose 
equation is given above. 

The fault-list computations for the eight possible patterns of Xo, X\, 
and X2 are summarized in Table XL 

As an example, consider the case *o = 0, X\ — 1, x 2 = l(line 4 in 
Table XI). We have 

x\ n = XoXi + X\X 2 + XoX 2 = a + b + c. 

The fault lists associated with the terms a = XqXi, b = X1X2, and c = 
xqx 2 for the specified values are: 

A=X nXi, B = X i UX 2 ; C = X nX 2 . 

Therefore, 

Xi„ = A n b n c = (Xo u Xi) n (Xi u x 2 ) n (X u x 2 ) 
= (X nx i )u (Xo n x 2 ) u (Xi n x 2 ) . 

4.2.4 Concurrent simulation 

In this case, each faulty signal value is computed separately. There- 
fore, high-frequency rejection can be performed on each signal individ- 
ually, using the equation for x in given in the preceding section. 

4.2.5 Suppression of short-duration detections 

We have considered the suppression of short pulses produced inde- 
pendently by each fault-free or faulty signal value. However, we did 
not consider the case of a short pulse of detection, when neither the 
faulty nor the fault-free signal incurs a pulse. This is illustrated by the 
case where Xo = 1, X\ = 1, x 2 = 0, for the fault-free signal and xo = 1, 
Xi = 0, x 2 = 0, for the faulty signal. This causes a short detection 
between U and t 2 . For deductive simulation, this short detection may 
be eliminated by using the formula 

Xm = (Xo n x,) u (Xo n x 2 ) u (X Y n x 2 ) 

independently of the fault-free signal pattern and after the high- 
Table XI — Fault-list equations for high-frequency rejection 

XO *l *1 *\n Xln 

(X n x,) u (X t n x 2 ) u (x„ n x 2 ) 
(Xo r\Xi) u iXi n x 2 ) u (Xo n x 2 ) 

(Xo n x^ u (X, n x 2 ) u (X n x 2 ) 
(X u n x t ) u (X, n x z ) u (X> n x 2 ) 
(Xo n x.) u (X, n x 2 ) u (Xo n X 2 ) 
(Xo n x,) u (X, n x 2 ) u (Xo n x.) 
(Xo n x.) u (X, n x 2 ) u (Xo n x 2 ) 
(Xo n x,) u (X, n x 2 ) u (X n x 2 ) 
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1 



frequency rejection has been performed. The term X is the set of 
faults detected at time U. 

The same method may be used for all the other simulation algo- 
rithms described earlier. 

4.3 Delay faults 

A fault that affects the transport delay associated with a signal is 
called a delay fault. Consider a fault that causes a delay to change 
from d to d'. When the signal at the site of such a fault changes, the 
signal value corresponding to the particular faulty circuit must be 
delayed by d' instead of d. 

4. 3. 1 Parallel simulation 

Two aspects of delay faults must be considered: injection of delay 
faults and the propagation of the effects of delay faults. Let us assume 
that a gate which is the site of a delay fault has been evaluated at time 
t, and the j th bit of the word represents the circuit with the delay 
fault. Let the normal and faulty delays be d and d', respectively, and 
let d < d'. Let x = (x u x 2 , ■ • • , x n ) and x' = (x'i, x'2, • • • , x' n ) be the 
vectors representing the old and new values, respectively. If for any 
i ¥> j, x'i ¥* Xi, a vector (x[, x'2, • ■ • , xj-i, xj, x'j+u • • • » **) will be 
scheduled to be applied to the gate output at time t + d. If x'j =£ Xj, the 
vector x' will be scheduled to be applied to the gate output at time 
t + d'. Similarly, if d' < d and x'j ¥= xj, the vector (xi, x 2 , • • • , Xj-i, x'j, 
Xj+i, • " ,x n ) will be scheduled for time t + d'. If x\ ¥* xi for any i ¥= j, 
then the vector x' will be scheduled to be applied at time t + d. The 
vectors for updating at the different times can be obtained from the 
old and new vectors by appropriate masks and logical operations. 

From the above discussion it should be clear that the effect of delay 
faults is to cause the signal values on the same lead in the presence of 
delay faults to change at different times. When one or more values in 
a vector change, the gates to which the signal fans out in the fault-free 
and all faulty circuits are scheduled for evaluation in parallel. There- 
fore, no special treatment is necessary for propagating delay faults. 
Since any signal change at the inputs of a device, faulty or fault-free, 
will cause an evaluation of the device (faulty and fault-free), delay 
faults will tend to increase the number of evaluations required. 

4.3.2 Three-list method 

The equations required for simulating delay faults using the three- 
list and deductive methods can be derived by treating the delay fault 
as an internal fault in a functional block. These equations can then be 
used for simulating delay faults without explicitly modeling them as 
faults in functions. 
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Let a be a delay fault which causes the delay associated with a signal 
to change from d to d'. Let f a be a fault variable, 19 which has the 
fault-free value of 0, but the fault a causes it to become 1. We shall 
represent the input and the output of the function used for modeling 
the delay fault by x and z, respectively. We assume that the evaluation 
is being done at time t = 0, and the value of the input x, t\ units of time 
earlier is represented by x(—ti). Two different cases must be consid- 
ered: 
Case 1: d < d'. The value of z at time d is given by 

if /„ then z = x(d — d') 
else z = x, 

which can be transformed into the equation 

2(d) = f n . X (d-d')+ f a .X. 

Case 2: d' < d. The value of z at time d' can be represented by a 
function as in Case 1, and transformed into the following equation: 

z(d') = f a -x+ f a .x(d'-d). 

The equations for the three-list method can be obtained from these 
equations using the method discussed in Section 2.2. 

Case 1: d < d' 

Z l (d) = [{a) H X\d - d')] U [X 1 f~! {«}] 

Z°(d) = [la) n A" ] U [X° n X°(d - d')] U [X°(d - d') n {«}] 



Z u (d) = [ZHd)UZ°(d)]. 
Case 2: d' < d 

Z\d') = [{a} n X 1 ] U \X\d' - d) n f«T] 

z°(d') = [T«T n.x°] u [x° n x°(d' - d)] u [XV - d) n {«}] 



Z"{d') = [Z'(d')UZV)]. 

4.3.3 Deductive simulation 

The functional equations derived in Section 4.3.2 can be used for 
deriving deductive fault lists for delay faults. The fault-list equations 
will depend on signal values as shown in Table XII. 

4.3.4 Concurrent simulation 

Since each fault is handled separately, the simulation of delay faults 
does not need any special processing. 
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Table XII — Fault-list equations for delay faults 



x x(d - d') z(d) 



Z(d) 



[X(d-d')n {a}]uLYn"pjj 
(«}]u|xn{^M 



[X(d - d') n {t 



xn 



R 



U [X(d - d') n {a}] U [XHX(d - d')] 



1 


1 


1 


[X(d-d')n{a}]U[XnX(d' 


-d)]U 


A"n"p| 




Case l:d<d' 


X 


x(d' - d) 


z(d') 


Z(d') 





[Xn{a)]U 



X(d' - d) n joj 



[xnx(d'-d)]u[xr\{a}]u 



X(d' -d)D 



H 



[X n {a}] u 



X(d' -d)n 



R] 



X(rf' - d) n |aj 



u[xn{a}]u[xnw-rf)] 



Case 2:d'<d 

4.4 Summary of results 

Changes in the fault-free and faulty values may occur at different 
times for the same line due to different rise and fall times and to delay 
faults. In the case of parallel simulation, a change in a single faulty or 
fault-free value on a line leads to computations involving the whole 
word (or pair of words). In the three-list and deductive methods, the 
addition or deletion of a single fault will require recomputation of 
complete lists. On the other hand, concurrent simulation treats each 
event, faulty or fault-free, independent of all other events and, there- 
fore, should require less computation. High-frequency rejection is also 
simpler in concurrent simulation than in the other methods. 

V. CONCLUSION 

We have compared parallel, multilist, deductive, and concurrent 
simulation methods with regard to their ability to simulate more than 
two logic values, different levels of simulation, and accurate timing 
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analysis. All the methods, except deductive, can handle any number of 
logic values without significant changes in the method. An extension 
of the deductive method to an arbitrary number of logic values is 
presented. Concurrent simulation appears to be the most convenient 
method of simulating an arbitrary number of logic values. 

All the methods, except concurrent, require the transformation of 
functional descriptions of high-level devices into Boolean equations. 
No such transformation is required for concurrent simulation. In fact, 
it is not even necessary to restrict operations in functional descriptions 
to Boolean operations if concurrent simulation is used. 

All the methods are capable of handling different rise and fall times, 
performing high-frequency rejection and simulating delay faults. Since 
concurrent simulation handles each event separately, these functions 
can be performed more easily and efficiently than the other methods. 

In addition to the aspects discussed here, two factors that must be 
considered in selecting a simulation method are storage requirements 
and speed. A detailed analysis of the speed and the storage require- 
ments of these methods is made in Ref. 11 based upon statistical data 
gathered from deductive simulation. 
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